Power fluctuations in a driven damped chaotic pendulum 
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In this paper we investigate the power fluctuations in a driven, dampted pendulum. When the 
motion of the pendulum is chaotic, the average power supplied by the driving force is equal to the 
average dissipated power only for an infinite long time period. We measure the fluctuations of the 
supplied power during a time equal to the period of the driving force. Negative power fluctuations 
occur and we estimate their probability. In a chaotic state the histogram of the power distribution 
is broad and continuous although bounded. For a value of the power not too close to the edge of 
the distribution the Fluctuation Theorem of Gallavotti and Cohen is approximately satisfied. 

PACS numbers: 05.45.Pq, 05.40-a, 05.45.Ac, 05.70.Ln 



Recently Gallavotti and Cohen (GC) Q derived a fluc- 
tuation theorem (FT) for chaotic systems. This theorem 
is concerned with the rate of entropy production a T av- 
eraged over an observation time r. It states that the 
ratio of probabilities of having a given entropy produc- 
tion a T to that of having (— a T ) is given by exp("rov). 
The theorem was proven for thermostatted Hamiltonian 
systems which are driven by external forces, under the 
condition that the system is sufficiently chaotic (of the 
Anosov type 0]) and the number of degrees of freedom 
is very large. All the systems considered had a certain 
form of time reversibility even though they were dissipa- 
tive. This was the result of the way the thermostat was 
imposed. The fluctuation theorem is actually a general- 
ization of the fluctuation dissipation theorem (FDT) to 
systems away from equilibrium i.e. driven systems in a 
steady state. Gallavotti || has shown that in the limit 
of vanishing driving force the FT reduces to the FDT. 

More recently Kurchan [Q and subsequently Lebowitz 
and Spohn || extended the FT to stochastic dynam- 
ics. Kurchan showed that the FT holds for finite sys- 
tems undergoing Langevin dynamics. The role of chaos 
is replaced by the stochastic fluctuations. Lebowitz and 
Spohn generalized the FT to general Markov processes. 
Particular role is played by systems displaying a form of 
local detailed balance. 

Our goal in this work was to test the FT for the case 
of a simple driven dissipative chaotic system that does 
not satisfy the conditions set by GC. For example, this 
system is most likely not of the Anosov type |J , and it is 
also not strictly thermostatted (energy is conserved only 
over a long period of time). Thus, it is deterministic but 
it is not time reversible because of the role of dissipa- 
tion. This example may serve as a representative that is 
more closely related to real life realizations. Surprisingly, 
(or not) we find that approximately, the FT is satisfied, 
taking into account that we have only three degrees of 
freedom, whereas the FT derived by GC should apply 
for the case of a large number of degrees of freedom. 

The system that we consider is that of a driven damped 



pendulum j7|-^0[. The driving force is periodic. Varying 
the strength of the driving force the pendulum can dis- 
play simple periodic behavior, a motion characterized by 
period doubling, tripling etc., or a fully chaotic motion. 
There is a lot of literature on such a system. What we had 
in mind was to consider the power bestowed on the sys- 
tem by the driving force during a complete period of the 
latter (or during an integer multiple of the period) . When 
the motion is periodic, the value of the power bestowed 
during any full period of the driving force is constant. 
But when the system is fully chaotic, there is a broad 
distribution and occasionally the power bestowed is neg- 
ative. We are going to investigate numerically the power 
distribution, and discuss the similarity to the GC Fluc- 
tuation Theorem. We are going to find that the power 
distribution is actually a good fingerprint to determine if 
the system is truly chaotic. 

We consider a planar pendulum subject to the force 
of gravity, a damping force resulting for instance from a 
viscous medium in which the pendulum is immersed, and 
a sinusoidal driving force acting around the pivot point. 
The equation of motion for the angle 9 is M: 
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Here b is the damping coefficient and Nd is the driving 
torque of angular frequency u>d- We are going to switch 
to dimensionless parameters. Denoting by uJq = g/i and 
to = 1/cl>0j we set t' = t/to and u = oj^/ujq. We also 
define 
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The equation of motion can now be written as a set of 
two first order differential equations (omitting the prime 
over the rescaled time variable): 



x(t) = y, 
y(t) = -cy 



sin(x) + Fcos(iut). 



(3) 
(4) 
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We have solved these equations numerically using 
Mathematica. We have used the values 



c = 0.05, oj = 0.7 



(5) 



and considered different values of the driving strength F 
in the range 0-1.0. Defining 



2tt 



Td 



(6) 



to be the period of the driving force, we will be interested 
in evaluating the quantity 



F 



t m +r 



dt y(t) cos(wt), 



with 
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(7) 



(8) 



and n is a fixed integer. We will be particularly interested 
in the case n = 1. (P) T is the average power supplied to 
the system by the driving force during n complete cycles 
of the latter. In the following we will calculate (P) Trn 
along a phase-space trajectory of the system, starting 
with m=100 to allow the transient to decay. We calculate 
consecutive averages over time periods of duration r and 
observe the fluctuations of the measured quantities. 
If we multiply Eq.(0) by y we get 
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Fy cos(wi) — cy 2 



Thus averaging over an arbitrary time period t 



-E) =F(ycos(ujt)) T 



(9) 



(10) 



where E is the total energy of the pendulum (kinetic plus 
gravitational). The first term on the right hand side of 
the equation represents the power input from the driving 
force and the second term represents the power output 
through dissipation with e.g. a viscous medium. The sec- 
ond term is always positive, whereas there is no restric- 
tion about the sign of the first term. We expect that in 
the limit of large r, the left hand side (= (E(t)-E(0))/t) 
will approach zero, if the system is in a steady state, and 
thus in that limit the two terms on the right should can- 
cel each other. For a simple harmonic oscillator, i.e. if 
the sinx in Eq.(Q) is replaced by x, then after the tran- 
sient term dies out the left hand side of Eq. ( |To| ) is zero 
for any r = n and 
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It has a maximum value of F 2 /(2c) when u> = 1 (at 
resonance) . This result is of course not applicable for the 
pendulum where the non- linear potential is present. We 
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FIG. 1. The average power vs. the magnitude of the driv- 
ing force. 



will see below that for the pendulum with fixed c and uj, 
(P) T behaves very differently as a function of F than this 
simple quadratic form. 

For small n, and in particular for n — 1, we will be 
interested in the fluctuations in (P) T , and in particular in 
its histogram. We will see that (P) T has a non- vanishing 
probability to be negative in the chaotic regime, and we 
will compare the probability of having (P) T = a, to the 
probability of (P) T = —a, within the valid range of a. 

We should also emphasize here the difference between 
our system and a thermostatted system considered in 
Ref. ( 0| ) . For a thermostatted system the dissipation is 
adjusted in such a way, by adjusting its velocity depen- 
dence, that the energy of the system remain identically 
constant, not only on the average over a long period of 
time, as is the case here. 

Varying the strength of the applied force the pendu- 
lum can display a simple periodic motion, a motion dis- 
playing periodic doubling, tripling, etc., or truly chaotic 
motion. In Fig. 1 we show the average power supplied 
by the driving force to the pendulum by using a time pe- 
riod of t = 400 Td for the averaging (after the transient 
component has practically died). The same results were 
obtained for the average power dissipated i.e. c(y 2 ) T . It 
is evident that beyond the initial segment F=0-0.4 which 
corresponds to periodic motion, the power curve is not 
smooth but vary significantly from point to point. The 
big jump in the power that occurs for F ~ 0.41 is associ- 
ated with the threshold for the pendulum to go over the 
top. The error bars are associated with chaotic motion 
when we do not have enough statistics to pinpoint the 
average power exactly. 

We now consider the histogram of the values of the sup- 
plied power during a time period r = t<j. When the mo- 
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FIG. 2. Histogram of the supplied power during a period Td 



FIG. 3. Plot of 7r(cr T )/7r(-<7 T ) for r = r d . The fit is to the 
form exp(Xr d a T ) as discussed in the text. 



tion is periodic, the histogram consists of a single peak. 
For example, for F=0.4 we get a single peak at =0.0134. 
For the case of F=0.5 there is period tripling, i.e. the 
Poincare plot shows three points. In that case, the his- 
togram consists of three peaks at P=-0.030, 0.128 and 
0.181. Thus the number of peaks in the power histogram 
coincides with the multiplicity of the period. For the case 
of fully chaotic motion the histogram becomes broad and 
continuous. We see that the histograms can be used as 
a good signature for chaotic motion. Only in the case of 
chaotic motion we obtain a broad, continuous distribu- 
tion of the power. Since we consider a system with a few 
degrees of freedom the distribution of power is bounded 
from below and from above. 

To plot the distribution of power we accumulated the 
data into bins. For F=0.6 we have chosen the bin size to 
be A = 0.036, and there are 21 bins spanning the range 
(-0.378,0.378). The bin counts for a run of 11300 periods 
of Td were found to be: 

{0, 0, 0, 0, 83, 305, 394, 677, 728, 648, 659, 795, 
1091, 1440, 1108, 1081, 729, 825, 429, 228, 80}. 

This is depicted in Figure 2. 

Using the histogram we calculate the ratio of 

ratio(i)= f ; e ^] + % i U.....6. (12) 
freq(ll - 1) 

The center of the ll'th bin corresponds to (P) T — 0. Let 
us define the effective temperature of the pendulum by 



-kT=-(y 2 ) =-(P) 
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(13) 



then define the rate of entropy production as 
{P) T _c{P) T 



kT 



(14) 



Since (P) T — A i, we can express ratio as a function 
of <7 T instead of i. In Fig. 3 we plot the probability 
of having an entropy production value a T during a time 
interval Td vs. entropy production — ay. We fit this to an 
exponential of the form 



exp(A T d oy), 



(15) 



where the average is over an infinite time period. We can 



with A a single fitting parameter. A = 1 provides a perfect 
fit for the first three points (solid line). A good fit to 
the first six points is obtained with A = 1.25 (dashed 
line). This is quite close to the value of A = 1 predicted 
by the FT for a system with a large number of degrees 
of freedom. Of course unlike the case for a system of 
many degrees of freedom the fit has to be cut off since 
the power distribution is bounded (the next data point is 
at infinity) . Thus if one is considering the properties of 
the histogram not too close to the edge, it approximately 
satisfies the FT. The deviation from A = 1 seems to be 
associated with the fact that the distribution is bounded 
due to the small number of degrees of freedom (finite size 
effect). 

We should mention that the dissipated power given 
by c(y 2 ) T also has a broad distribution for the case of 
chaotic motion, but of course it is never negative. Thus 
unlike the thermostatted case the two histograms of the 
supplied power and dissipated power are different. 

We also calculated the histogram of supplied power 
for a time period r = 2Td, see Fig. 4.. The sum of 
frequencies is 10,000 and the bin size is 0.027. In that 
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FIG. 4. Histogram of the supplied power during a period 
2r d 



case the distribution is narrower as expected, although 
there are still periods of negative power. Of course in 
the limit of large r the histogram will be peaked at the 
average value. 

It is now clear that chaos is characterized by a wide, 
continuous distribution of entropy production rates, 
whereas in the non chaotic regime the histogram is char- 
acterized by discrete peaks whose number is equal to the 
multiplicity of the period in terms of the period of the 
driving force. Thus the power distribution is a good sig- 
nature for the chaotic state. The supplied power can be 
negative for time periods equal to the period of the driv- 
ing force, but the the distribution is centered about a 
positive average value. For values of the power not too 
close to the edge the ratio of probabilities 
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This suggests that at least near the center one may ap- 
proximate the power distribution by a gaussian of the 
form 



7r T (P) oc exp 



4(p>: 



(P-(P)oo) 



(17) 



It should be emphasized that a gaussian is only one pos- 
sible solution satisfying the FT, and other solutions de- 
viating from a gaussian are possible. 

If we consider N identical (uncoupled) driven dissipa- 
tive pendulum, then the total power supplied which is 
the sum of individual powers will have a gaussian distri- 
bution in the limit N — * oo as predicted by the central 
limit theorem. One can consider N coupled pendulums, 
and in that case it will be very interesting to check for 



the validity of the fluctuation theorem when the motion 
of the system is chaotic and N — > oo. In that case the 
power distribution is expected to be unbounded. It need 
not necessarily be gaussian. Devaitions from gaussian 
behavior were observed for turbulant flows [Q. 

Gallavotti and Cohen used the time reversal invari- 
ance of the particular thermostatted systems that they 
considered in order to prove the FT. Such a symmetry is 
not present in our system, although energy is still con- 
served when averaged over an infinitely long time period. 
Of course it will be very interesting to find a proof of 
the theorem for a steady state of a dissipative chaotic 
system for which the energy is conserved only over long 
time periods if indeed the theorem is still true. 

Another ingredient of the proof of GC is the validity of 
the SRB measure on the attractor. This measure implies 
that the relative time the system spends at a certain re- 
gion of the attractor is inversely proportional to the pos- 
itive expansion rate at that region as given by the local 
positive Lyapunov exponents (for a long trajectory cen- 
tered at that local). We carried out jllj a preliminary 
investigation of the attractor of the chaotic pendulum 
and showed that there is a positive correlation (although 
not perfect) between the time spent by the system at a 
region of the attractor (divided into cells) and the in- 
verse local expansion rate as given by the local positive 
Lyapunov exponent. 
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